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ABSTRACT 


The mutual impedance of nonplanar-skew sinusoidal dipoles is 
presented rigorously as a summation of several exponential integrals 
with complex arguments. 
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I . INTRODUCTION 

In 1957, H. E. King [l] published the mutual impedance 
expressions for parallel dipoles in terms of sine-integrals and 
cosine-integrals. This paper gives the analogous expressions for 
non-parallel dipoles. 


Previous communications consider the mutual impedance of 


coplanar^s 
dipoles [3 


cew sinusoidal dipoles 
and coplanar V dipoles 


planar and nonpl anar-skew 
A computer program is 


available [5] for thin-wire antennas" and scatterers with sinusoidal 
bases and Galerkin's method. In this program, subroutine GGMM 
evaluates the mutual impedance between filamentary monopoles with 
sinusoidal current distributions. The mutual impedance of two V 
dipoles is the sum of four monopole-monopole impedances. This paper 
presents the expressions programmed in GGMM. For thin-wire structures, 
the speed and accuracy [6] of the sinusoidal-Galerkin program is 
attributed largely to the use of these expressions. 


When the monopoles are relatively far apart, the impedance can 
be calculated faster with numerical integration (as in subroutine 
GGS) . The expressions given here are most useful when the monopoles 
are close together. 


II. NEAR-ZONE FIELD OF THE SINUSOIDAL MONOPOLE 

Many books [7] and papers [8] give the near-zone fields of a 
sinusoidal line source in free space. This section presents the 
fields of a time-harmonic electric line source (monopole) in a homo- 
geneous conducting medium. The filamentary source is located on the 
z axis with endpoints z x and z 2 and length d as in Figure 1. A con- 
venient expression for the current distribution is 

I, sinh(z 2 - z) + I 2 sinh(z - z 2 ) 

(1) 1(2) = 

sinh d 


(2) y = s /ye" 


where the complex constants Ij and I 2 denote the endpoint currents. 
The complex frequency is s = jw, and w may be real or complex. In 
Eq. (1) and the subsequent equations, quantities of the form yx are 
denoted simply by x. This convention applies only when x represents 
a linear dimension (such as d) or a metric coordinate such as z. 
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Fig. 1--Electric monopole and the coordinate system. 


When the sinusoidal monopole radiates in a homogeneous conducting 
medium with complex permeability u and complex permittivity e, the 
field is given by 


(3) 


„a r 

= n- i 

p L 


-Rn 


-R, 


(4) 


E z = nA 


(5) H 


( I i e - I 2 e ) sinh d + (I x cosh d - I 2 ) e cos ^ 

-R 2 

+ (I 2 cosh d - I x ) e cose. 


- R2 -Ri 

(Ij - I 2 cosh d) — + (I 2 - I x cosh d) 

r 2 Ri J 


-Ri 


= |^ ( I sinh d cose 1 + I x cosh d - I 2 ) e 

-R 2 -| 

- (I 2 sinh d cos ^ - I 2 cosh d + I x ) e 
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^ A (4tt slnh d) 


(7) n = /y/e 


where (p, <^z) denote the cylindrical coordinates. (R . , e. , <j») are the 
coordinates in a spherical system with origin at the endpoint . 

The sinusoidal monopole has point charges at the endpoints, but they 
will disappear when another monopole is connected to form a dipole. 
For this reason, the contributions from the point charges are omitted 
in the above field expressions. 


III. MUTUAL IMPEDANCE OF SINUSOIDAL MONOPOLES 

In the induced emf formulation, the mutual impedance of coupled 
dipoles is 


(8) Z = - L,{t) • Ej(t) dt 


where I 2 (t) denotes the current distribution (normalized to unit 
terminal current) on dipole 2, and E^t) is the field of dipole 1 when 
it transmits with unit terminal current. Distance along the axis of 
dipole 2 is denoted by the coordinate t. In Eq. (8), Ej may be ex- 
pressed as the sum of the fields from each of the monopoles comprising 
dipole 1. Furthermore, the integral in Eq. (8) is the sum of the 
integrations over each of the monopoles comprising dipole 2. Thus, 
the dipole-dipole mutual impedance may be expressed as the sum of four 
monopole -monopole impedances. 

With no loss of generality, we again locate the source monopole 
on the z axis with endpoints at Zj and z 2 and length dj = z 2 - Zj. As 
shown in Figure 2, the receiving monopole is located on the t axis 
with endpoints at and t 2 and length (fc = t 2 - tx . Monopole 2 lies 
in the plane y = d. The coordinate origin for (x,y,z) is located at 
the apparent intersection of the t axis and the z axis in Figure 2. 

The origin for the t coordinate is at (x»y,z) = (0,d,0). Monopole 1 
is a segment of dipole 1, and monopole 2 is a segment of dipole 2. 

The monopole -monopole impedance will not necessarily satisfy the 
reciprocity theorem. The reference directions for current are z and 
t. The mutual impedance is defined by 

(9) Z - -J^ 2 T(t) E(t) dt 
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Fig. 2— Monopole 1 lies on the z axis and monopole 2 on the 
which is parallel with the xz plane and displaced a 
distance d from it. 


(10) E(t) = E p cos <f) si nip + E z cosiJj 

(11) cos $ = x/p 


In Eqs. (8) and (9) the integration variable is 
current distribution on monopole 1 be given by Eq 


t and not yt. 
and 


( i \ 


j -u 
wi< 


on monopole 2 is 

(12) T(t) = [ Tj sinh(t 2 - t) + J 2 sinh(t - t ^ /sinh d 2 
where and \ 2 are the endpoint currents. 


t axis 


Let the 
6 current 
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To be consistent with Eqs. (8) and (9), the currents I{z) and 
I (t ) must have unit value at one end of the monopole and vanish at the 
other end. Thus a double -subscript notation Z-j j is convenient for the 
mutual impedance of two monopoles. The first subscript i has value 
1 if Ij = 1 and value 2 if I _2 = 1. The second subscript j has value 
1 if Ii = l‘and value 2 if I 2 = 1. (The monopole has unit current at 
the endpoint corresponding to the dipole terminals.) The first sub- 
script is associated with the "source monopole" whose field E(t) 
appears in Eq. (9). The second subscript is associated with the 
"receiving monopole" with current I (t). As illustrated in Figure 3, 
the source monopole is always a segment of the source dipole. 



DIPOLE DIPOLE 


Fig. 3--The mutual impedance of coupled dipoles may be expressed 
as the sum of four monopole-monopole impedances: 

Z= z n + Z 12 + Z 21 + Z 22 . 


Omitting the lengthy derivations, the mutual impedance Z-j j for 
skew monopoles is J 


( 13 ) 










+ 


z 


e 


m 



5 



(14) B = 


n 


(16 it sinh d 1 sinh d 2 ) 

where m = 2/i and n = 2/j. The functions are defined by: 

AZ-j COSiJi 

(15) F ik = 2 sinh d l e E(R 1 * + £z.j cos^ - At) 

where 

. . 1 / 

A = (-1) . The functions G. are defined as follows: 

k£ 

(16) = E(R 2 + mz 2 + nt - jg) + E(R 2 + mz 2 + nt + jB) 

-EfRj + mz 1 + nt - jB) - E(R X + mz x + nt + jg) 

where m = (-1 ) k , n = (-1) £ and p = nt + nc. 

In Eq. (16) and all the subsequent equations, j = v^T. Let 

(17) b = c cost 

(18) c= d/sint 


In Eqs. (15) and (16), t denotes the position of an observation point 
somewhere on monopole 2.< R x and R 2 are distances from the endpoints 
of monopole 1 to this observation point. Finally, the E functions are 
defined as follows: 


(19) 


E<« + 

TiD 


yw 


dw 


. D 
l ' O M 


W 


where a and B are real quantities with dimensions of length, a is a 
function of t, = a(ti), a 2 = a(t 2 ) and y is given by Eq. (2). In 
the right-hand side of Eq. (19), we abandon the convention of letting 
x represent yx. 


The integral in Eq. (19) may be expressed in terms of exponen- 
tial integrals as follows: 
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= Ei(vj) - E 2 ( v 2 ) + j 2 n-rt 


(20) s = S —~ 

The integration path is a horizontal line in the w plane or an inclined 
straight line from v^ to v 2 in the v plane. The integer n is zero 
unless this path intersects the negative real v axis at a point between 
vi and v 2 . When there is such an intersection, 

a) n = 1 if Vj lies above the real axis and v 2 is below, or 

b) n * -1 if vj lies below the real axis and v 2 is above. 

A subroutine EXPJ is available for the integral in Eq. (19). This 
subroutine is used with GGMM for skew monopoles and dipoles. Although 
the expressions given here fail for parallel monopoles, the subroutines 
do not. 


IV. CONCLUSIONS 

Rigorous expressions are presented for the mutual impedance of 
nonpl anar-skew sinusoidal monopoles in a homogeneous conducting 
medium in the complex frequency domain. 


7 



APPENDIX 


Two basic integrals are required to derive the results presented 
herein. These "reaction integrals" were evaluated by N.H. Geary in 
October 1968 and April 1969 with substitution of variables and partial 
fraction expansions. The integrals are defined as follows: 

(22) P = 2 J^ 2 


(23) Q = 2 sinifi J^ 2 e~ YR dt 

In this appendix, a subscript i is understood on the quantities R, 

6 and z. The subscript may take on values 1 and 2. In the integrands, 
all quantities are functions of t except Y and 'I'. In terms of the 
E function in Eq. (19), P and Q are given rigorously as follows: 


_±Yt -YRr COSif; 

L R 


cos 4 1 cos 6 s 


m 


dt 


(24) 

P = 

-YZ 

e 

E(Xj) 

N 

>- 

0) 

1 

E(x 2 ) 

+ e' YZ E(x 3 ) 

- e YZ E(x 4 ) 

(25) 

Q = 

-YZ 

e 

E(Xi) 

+ e YZ 

CM 

X 

LlJ 

+ e' YZ E (Xj ) 

+ e YZ E(x 4 ) 



- 2 « 

.*« c 03 * E(x 

5 ) 



(26) 

x i : 

» R + 

t - z 

+ jb ± 

jc 



(27) 

X 2 : 

= R + 

t + z 

- jb + 

jc 



(28) 

*3 : 

= R + 

t - z 

- jb + 

jc 



(29) 

x 4 5 

■ R + 

t + z 

+ jb + 

jc 



(30) 

X 5 ! 

= R + 

t ± 2 

COStJ> 





All quantities not defined here are defined in the main text. 
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